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Deep  Level  Defects  in  Semiconductors 


Abstract 

Using  Racah's  irreducible  tensor  operator  formalism,  a  generalized  and 

more  refined  treatment  of  the  d-electron  matrices  for  transition  metal 

defects  in  crystal  has  been  given.  The  resulting  Coulomb  and  exchange 

interaction  matrices  have  been  used  to  calculate  the  electronic  structures 

of  6aAs:Cr  and  GaAs:Cr  ancf  interpret  the  optical  data  on  MgF2:Co  and 
2+ 

MgF2:Mn  .  The  significance  of  the  new  theory  is  explained.  From  the 

photoluminescence  and  optical  absorption  data,  the  crystal  field  parameters 

have  been  derived  for  GaAs:Fe  and  InP:Fe.  The  theory  of  the  charge  transfer 

2+ 

effects  on  the  electronic  structure  of  GaAsiCr  have  been  discussed  and  a 

ligand-field  model  has  been  developed  for  explaining  the  optical  and  EPR 
2+ 

results  of  GaAsrCr  . 


1.  Introduction 


The  program  described  in  this  final  report  is  part  of  an  over-all 

program  dealing  with  defect  centers  in  III-V  semiconductor  materials 

that  are  of  importance  for  opto-electronic  device  applications. 

Specifically,  the  investigations  relate  to  the  optical  studies  of 

transition  metal  impurities  in  GaAs  and  InP  giving  rise  to  deep  level 

centers.  The  report  consists  of  the  results  on  (i)  the  derivation  of  the 

crystal  field  parameters  that  will  best  explain  the  photoluminescence 

results  of  InP:Fe  and  optical  absorption  at  low  temperatures  of  GaAs:Fe; 

(ii)  the  derivation  of  the  generalized  d-electron  interaction  matrices 

for  defect  centers  in  crystals  and  the  calculations  of  the  electronic 

structures  of  GaAs:Cr  and  GaAs:Cr  ;  and  (iii)  the  considerations  of  the 

charge  transfer  effects  affecting  the  g-factors  and  the  zero  field 

2+ 

splitting  parameters  in  GaAs:Cr  leading  to  the  interpretation  of  the  EPR 

results  of  the  system.  In  addition,  from  the  generalized  treatment  of  the 

electron-electron  interaction  matrices  we  have  also  provided  the 

2+  ?+ 

interpretation  of  optical  data  on  MgF2:Co  and  MgF^rMn  .  The  remaining 

part  of  this  report  consists  of  four  parts  (Section  2-5)  and  supplements 

1  2 

our  published  article,  ’  as  well  as  the  presentations  at  the  American 
Physical  Society  meetings  and  the  ones  that  have  been  submitted  to 
journals  or  are  under  preparation.  While  this  report  summarizes  the  results 
of  our  studies  thus  far,  further  information  on  any  of  the  topics  discussed 
in  this  report  may  be  obtained  from  the  investigators. 


2.  Intracenter  Transitions  in  InP:Fe  and  GaAs:Fe 

The  transition  metal  impurities  associated  with  compensation  and 

control  of  electrical  properties  of  semiconductors  give  rise  to  deep  levels 

in  the  forbiddden  gap.  While  the  nature  of  the  deep  level  impurities  in 

II-VI  semiconductors  has  been  studied  extensively  by  contrast  the  nature  of 

these  impurities  in  III-V  compounds  has  not  received  much  attention  until 

recently.  The  positions  of  the  impurity  centers,  their  charge  states,  and 

the  nature  of  the  recombination  processes  are  all  very  relevant  to  the 

behavior  of  III-V  materials  and  their  uses  in  device  applications. 

The  present  study  on  the  interpretation  of  the  photoluminescence  and 

optical  absorption  experimental  results  was  undertaken  as  part  of  our 

systematic  study  of  the  deep  level  defect  centers  in  GaAs  and  InP.  The 

studies  on  Fe  in  III-V  materials  have  been  limited  to  GaAs  in  most 
7-12 

studies  .  The  only  studies  of  Fe  in  InP  are  the  recent  measurements  of 
12  13  2+ 

Koschei  et  ^  ’  establishing  Fe  as  a  possible  charge  state  for  the  ion 
The  photoluminescence  emission  peaks  observed  in  the  infrared  and  their 
behavior  with  temperature  have  been  used  in  our  investigation  to  establish 
that  the  mechanism  can  be  explained  as  arising  from  intracenter  transtiton. 
There  are  four  prominent  lines  observed  in  the  photoluminescence  spectrum 
of  InP:Fe.  They  are  at  2801  cm~^,  2819  cm~^,  2830  cm'^,  and  2843  cm~^.  A 
comparison  of  the  PL  spectrum  with  the  optical  absorption  indicates  that 
crystal  field  theory  applied  to  Fe  (3d  )  in  a  tetrahedral  (T^)  environment 
of  the  host  crystal  InP  can  satisfactorily  account  for  the  observed  lines. 
Group-theoretical  arguments  lead  to  the  splitting  of  the  term  D  of  the 
free  ion  Fe^^  into  orbital  doublet  ^E(r3)  4nd  orbital  triplet  ^^2^ ‘6'  ^d 
Due  to  the  second  order  spin-orbit  interaction  effects  the  orbital  doublet 
will  split  into  five  levels  and  the  orbital  triplet  will  give  rise  to  six 


levels  according  to 


Dj  X  Fj  =  Fj  +  rj  t  r,  +  r4  +  Fj 

and 

DjX  s'l-l^rj^Fj  +  F^  +  Fj 

(see  scheme  in  Fig.  1).  Due  to  lower  nonequilibrium  electron  populations 

5 

of  the  higher  levels  of  the  excited  state  J2  and  the  possible  rapid 

5 

thermal izati on  of  higher  states  to  the  lowest  level  of  by  non- 

radiative  transitions,  one  need  consider  only  the  transitions  from  the 
5  5 

lowest  level  of  I2  to  the  split  levels  of  E.  The  transition 
fc(To)  -»■  ro(E)  is  forbidden  and  this  leaves  four  allowed  transitions. 

be  c 

2+ 

In  Table  I  we  have  summarized  the  frequencies  observed  for  InP:Fe  , 
the  assignments,  and  the  expressions  for  the  energy  differences  in  terms 
of  the  crystal  field  parameters  (=10Dq),  Xj  and  X^  (the  spin-orbit 
interaction  parameters).  From  the  expressions  for  the  energies,  the 

separation  between  the  observed  PL  peaks  should  be  equal  to  the  theoretical 

2  -12 
splitting  6X2/ A  .  Using  an  average  value  of  14  cm  for  this  separation  bX^/A 

from  the  experimental  data,  a  value  of  84  cm"^  is  derived  for  x^ 

corresponding  to  300  cm*'  Dq.  For  the  range  of  260  to  310  cm'^  for 

Dq,  one  gets  X2  =  81±4  cm”  .  The  mean  values  of  the  crystal  field 

parameters  derived  in  our  study  are 

A  =  3028  cm“^  and  Jx^l  =  X2  =  83  cm”^, 

g 

Table  II  shows  the  data  on  low  temperature  optical  absorption  by 

GaAs:Fe.  The  assignments  of  the  observed  transitions  by  Ippolitova 
g 

,  are  also  shown  in  Table  II.  The  labels  correspond  to  the  numbers 

given  in  Fig.  1.  On  the  basis  of  the  assignments,  the  following  crystal 

g 

field  parameters  have  been  obtained  by  the  previous  authors. 


I 

Table  I 


fa 


Intracenter  Transitions  of  InP:Fe 


2+ 


Transition* 


Frequency ( cm” 


Energy 


*  ' 

^  • 

Transition* 

Frequency(cm~^) 

Revised 

Energy 

;  (Previous  Assignment) 

Assignment 

r 

m  * 

1  -  6 

3002 

1  -  6 

A 

+ 

3Xj^  138  X^/Sa 

’* 

2 

vV 

1 

1  -  9 

3092 

2  -  8 

A 

- 

2X  +  132  ^2 

5  A 

\r' 

2-6 

2988 

2  -  6 

A 

+ 

3X2  +  108  X^/Sa 

i 

,2 

* 

3044 

4  -  8 

A 

+ 

^2 

•  • 

^  2 

2  -  8  J 

, 

30X^ 

i 

A 

+ 

h"  A 

«  , 

3  -  6 

2979 

3  -  6 

A 

+ 

3X^  +  78  X^/Sa 

3  -  6 


2979 


3  -  6 


o 


A  =  2995±10  cm"^,  =  18±2  cm~^,  =  77±6  cm~^ 

We  have  made  a  systematic  review  of  the  above  optical  absorption  data  and 

the  results  for  the  crystal  field  parameters  obtained  by  the  previous 
g 

authors.  Though  the  covalency  effects  can  reduce  the  free  ion  values  for 

14 

and  X2  by  'v-20%  as  pointed  out  by  Wong  it  seems  highly  improbable  that 
a  value  as  low  as  18  cm"^  can  be  obtained  for  x^.  There  seems  to  be  an 

9 

error  also  in  the  sign  of  X^.  These  prompted  our  review  of  the  assignments 
of  the  observations,  particularly  those  at  3092  cm  and  3044  cm”^.  The 
assignments  of  the  transitions  to  1-9  and  2-7  or  2  .especti vely  do  not 

seem  plausible  on  the  basis  of  the  value  of  Xj^  or  ^parations  of  the 
observed  transitions.  A  revised  assignment  of  these  two  bands  to  2-8  and 
4-8  transitions  seem  to  be  reasonable  and  these  are  allowed  by  selection 
rules.  With  these  revised  assignments,  all  the  observed  absorption  lines 
for  GaAs.'Fe  yield  the  following  values: 

A  =  3028  cm'^,  X^  =  -44  cm”^  and  X2  =  B3  cm’^. 

We  have  compared  the  values  for  InP:Fe  and  GaAsrFe  in  Table  III,  and  there 
is  a  close  agreement  between  the  two  sets  of  values  and  the  values  obtained 
by  Bykovski  e^al^.  who  report  a  value  of  -90  cm”^  for  X^.  Thus  in 
summary,  the  photoluminescence  data  on  InP:Fe  and  the  low  temperature 
optical  absorption  data  on  GaAs:Fe  have  been  interpreted  as  arising  from 
intracenter  transitions  and  the  crystal  field  parameters  of  A,  x^^,  and  X2 
have  been  derived  from  the  above  data.  The  revised  assignments  for  GaAs:Fe 
provided  by  our  study  give  better  sets  of  values  as  may  be  seen  from  Table 
III  and  by  comparison  with  results  on  II-VI  materials. 


System 


A  =  10  Oq 


InP:Fe^^ 

3028 

GaAsiFe^^ 

3028 

(2995) 

Values  in  Parentheses  are  from  Ref.  9 


3.  Energy  Levels  of  GaAstCr^^  and  GaAs:Cr^^ 

In  this  section  we  present  a  major  advance  for  the  first  time  in  the 
theoretical  treatment  of  the  d-electron  interactions  in  order  to  facilitate  the 
theoretical  understanding  of  the  experimental  observations  on  the  III-V 
materials  of  GaAs  and  InP  containing  transition  metal  impurities  as  defect 
centers.  It  is  well  known  that  up  to  now  most  of  the  theoretical 
treatments  have  been  restricted  to  the  use  of  the  "pure"  d-orbitals  in 
calculating  the  electron-electron  interactions  in  complexes  and  solids.  In 
our  attempts  to  investigate  systematically  the  nature  of  defect  centers  in 
GaAs  and  InP  we  have  encountered  the  difficulties  associated  with  the 
unavailability  of  the  explicit  generalized  form  of  Tanabe-Sugano  matrix 
elements  for  the  single  particle  electronic  wavefunctions  which  are  not 
of  the  "pure"  d-character.  The  theory  that  was  widely  adopted  for 
interpreting  the  optical  spectra  of  transition  metal  ions  in  solids  and 

It 

complexes  is  the  one  by  Racah  using  only  three  parameters  B,  C,  and  a. 

As  is  well  known,  despite  the  useful  simplifications  this  theory  provided, 
it  suffers  from  serious  drawbacks  in  neglecting  the  important  solid  state 
effects  such  as  the  modifications  of  the  wavefunctions  from  the  pure  d 
character.  There  are  a  number  of  examples  in  which  it  has  been  shown  that 
this  simple  theory  is  glaringly  inadequate.  Therefore,  it  is  extremely 
important  to  improve  the  theory  and  we  have  adopted  the  strong  field 
coupling  scheme  ’  (since  many  physical  cases  of  interest  come  under  this 
category),  to  provide  the  refinement  to  d-electron  interactions. 

In  our  treatment,  which  also  adopts  Racah's  irreducible  tensor  operator 
formalism,  we  have  derived  the  generalized  electron-electron  interaction 
matrices  in  terms  of  ten  independent  Coulomb  and  Exchange  parameters  Vj^,  V^, 
...V^Q  (or  a,b,c,...j)  defined  in  the  standard  Dirac  notation  and  the 
one-electron  orbitals  ",  -  ,  ri,  and  r.  in  the  cubic-field  representations. 


2  3  4  5 

The  matrix  elements  for  d  ,  d  ,  d  ,  and  d  configurations  in  octahedral 

synr.ietry  are  tabulated  in  the  attached  reprint  (Ref.  2).  As  might  be 

inferred,  if  i  ,  ,  t;,  are  reduced  to  pure  d-type  orbitals,  these 

15 

matrix  elements  reduce  to  the  forms  given  by  Tanabe  and  Sugano.  In 
addition,  using  group-theoreti cal  methods  the  matrix  elements  of  the 
complementary  states  (n  =  2,  3,  4,  5)  have  been  derived  by  us.  Of 

particular  significance  in  the  general  expression  that  we  have  derived 
connecting  the  complementary  states  in  the  form 

<'t^‘"'(S^r^)e'^'^  (52r2)  Ke^/ r^2) 

e''■^(S2^2)SlM^^^,)>  -  <t^'(S^r^)e^(S2r2) 

SrM^Mj,  l(e2/r^2)  I  ^2  {S2r2)SrM^Mj,) 

=  (3-m)V^  +  4(3-m)V2  +  (24-4m-6Ji)  (V3  +  V^)  + 

+  (5t-10)Vg  +  (2m  +  3i-  -  12)V7  +  (2m  +  3l  -  12)Vg/v^  -  2(3-m)VjQ 

where  the  integers  m  and  £  are  connected  to  n  (in  d^  configuration)  by 
the  relation 

m  +  i  =  n  or  10-n. 

Equally  significant  is  the  fact  that  we  have  shown  that  the  diagonal  matrix 

elements  of  states  differ  from  the  corresponding  elements  of  d*^  by 

different  amounts  that  are  linear  combinations  of  V/s.  This  is  in  contrast 

to  the  constant  differences  resulting  in  the  case  of  pure  d-orbitals. 

3  4 

The  matrix  elements  derived  in  the  present  study  for  d  and  d 

3+  2+ 

configurations  were  used  to  obtain  the  energies  for  GaAs:Cr  and  GaAs:Cr  . 
Recently,  Hemstreet  and  Dimmock^^  have  calculated  the  above  electronic 
structure  by  modifying  the  e  and  t2  symmetry  orbitals  using  parameters  R^^, 


12 


Rtt>  Rgt  deduced  from  calculations.  We  have  calculated  the  energies 
2+  3+ 

of  GaAs:Cr  and  GaAs:Cr  using  our  generalized  matrix  elements  and  these 

results  are  presented  in  Figures  2  and  3.  Except  for  some  minor  differences 

2 

that  we  have  pointed  out  in  our  publication  ,  there  is  general  agreement 

18 

between  our  results  and  those  of  the  earlier  study. 

While  our  new  theory  has  more  parameters  (eleven  of  them  Vj^,  ^2  ••• 

VjQ,  a)  instead  of  four  parameters  (A,  B,  C,  a)  or  more  precisely  only 
three  parameters  (B,  C,  a)  used  by  Racah^^,  the  full  significance  of  our 
theory  lies  in  the  potential  application  to  a  variety  of  physical  problems 
in  the  correct  way  and  the  generality  of  our  superior  formalism.  This  is 
further  clarified  in  Part  5  of  this  report  involving  refined  interpretation 
of  the  spectra  of  other  systems  and  the  evaluation  of  Coulomb  and  Exchange 
interaction  constants. 


Energy  Levels  for  GaAsiCr 
(a)  Previous  study  (Ref.  18)  (b)  Present  study 
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2+ 

4.  Theory  of  g-factors  and  Charge  Transfer  in  GaAs:Cr 

As  part  of  our  continuing  studies  of  transition  metal  impurity  defect 

centers  in  GaAs,  in  this  section  we  present  the  results  of  some  of  our 

2+ 

preliminary  investigations  on  the  nature  of  Cr  center  in  GaAs  and  the 

2+ 

interpretation  of  EPR  experimental  data.  Although  GaAs:Cr  has  been  the 

subject  of  several  theoretical  and  experimental  studies  particularly  with 

respect  to  the  electronic  structure  of  the  ground  state  of  the  impurity, 

many  theoretical  questions  for  the  excited  states  of  the  impurity  ion 

19 

still  remain  unanswered.  Electron  paramagnetic  resonance  (EPR)  and 

20  2+ 
ultrasonic  attenuation  experiments  have  established  that  GaAs:Cr 

undergoes  a  Jahn-Teller  distortion  resulting  in  a  change  from  tetrahedral 

(Tj)  to  tetragonal  (02^)  symmetry  at  the  Cr  site.  The  upper  state  E 

5  5 

of  the  0.84  eV  transition  (  £-»•  T2)  observed  in  photoluminescence 
21  5  5 

experiments  splits  into  A^  and  states  of  D2J  symmetry  while  the 

ground  state  T2  splits  into  82  and  ^E  states  in  the  tetragonal  symmetry 

as  a  result  of  Jahn-Teller  (JT)  distortion.  (See  Fig.  4)  Such  results  are  also 

22  23 

consistent  with  photconducti vity  measurements.  ’  In  a  similar  study  on 

2+  24 

analogous  II-VI  materials  with  Cr  ,  Vallin  and  Watkins  approached  the 

problem  through  a  molecular  orbital  treatment.  Though  they  used  only  one 

parameter,  the  study  did  reveal  the  importance  of  the  charge  transfer 

effects  on  the  magnetic  properties  of  the  impurities  in  such  systems. 

In  view  of  the  presence  of  Jahn-Teller  distortion  and  the  importance 

of  charge  transfer  effects  we  have  studied  as  part  of  this  project  the 

2+ 

molecular  cluster  Cr  -As^.  The  present  molecular  orbital  treatment  in  the 

framework  of  the  ligand-field  theory  involves  the  construction  of  molecular 

2+ 

orbitals  using  the  d-orbitals  of  Cr  and  the  outermost  s  and  p-orbitals  of 
the  As  atoms.  In  terms  of  the  cubic  field  representation  0,  i ,  n,  ’ . 


i; 


cluster  wavefunctions  v  ,  i'  belonging  to  and  ’i.,;  ,i.  belonging  to 
and  were  constructed  and  the  expressions  for  g^^,  gj^,  and  D  in  the  spin- 
Hamiltonian  were  derived  in  our  study.  The  details  of  such  calculations  and 
the  exact  forms  of  the  expressions  for  wavefunctions  involving  the  admixture 

parameters  will  be  given  in  an  article  to  be  published.  The  EPR  results 

19  -1 

obtained  by  Krebs  and  Stauss  give  g^,  =  1.974,  g^^  =  1.997,  and  D  =  -1.860  cm  , 

2+ 

In  our  study  we  have  obtained  the  spin-orbit  matrix  elements  ^  (Cr  )  and 

Cp  p  303  cm'^  and  1246  cm"^  respectively.  These  values  lie  close  to 

the  corresponding  values  of  236  cnf^  and  1273  cnf^  given  by  Abragam  and 
25 

Bleaney  .  In  the  present  investigation,  the  overlap  and  two-center  integrals 
were  calculated  and  the  values  of  admixture  parameters  (x's)  that  will  best 
explain  the  EPR  results  (gjp  gj^,  and  D)  were  derived.  These  values  are 

^dp„  =  0-^91 

Xds"  =  0.281 

=  0.110. 

dp„ 

Assuming  the  X's  are  sums  of  the  corresponding  overlap  and  charge  transfer 

covalency  parameters,  the  values  of  the  charge  transfer  covalencies  obtained 

in  the  present  study  are 

r.„  =  0.583 
dp 

rds"  =  0.214 
r.„  =  0.049. 

71  2+ 

These  values  are  significant  and  reveal  the  nature  of  the  defect  center  Cr 

in  GaAs  and  can  also  be  used  to  derive  more  accurate  Coulomb  and  exchange 

2+ 

matrix  elements.  The  cluster  Cr  -  As^  is  shown  in  Fig.  5. 


5.  Interaction  Constants  For  Ions  In  Crystals 

As  was  mentioned  earlier,  one  of  our  main  goals  was  to  provide  a  refined 

treatment  for  the  interaction  matrices  for  transition  metal  ions  in  crystals. 

In  the  process  of  deriving  a  new  formalism  and  providing  the  generalized 

matrices  we  have  come  across  the  detailed  optical  spectra  of  systems  involving 

2+  2+ 

transition  metal  ions  Co  and  Mn  in  MgF2  system.  Though  the  system  is  not 

of  direct  interest  to  our  series  it  was  felt  that  the  new  theory  should  be 

verified  by  application  to  a  system  where  there  are  abundant  lines  observed 

experimentally.  The  experiments  on  the  above  mentioned  system  have  made 

available  many  sharp  lines  of  the  transition  metal  ions  in  crystals  by  lifting 

the  spin-forbidden  characteristics  of  the  excitations  through  exchange 

interactions  with  color  centers.  The  refined  formalism  that  we  have  derived 

2+ 

was  applied  in  the  present  investigation  to  analyze  the  spectra  of  MgF2:Co 
and  MgF2:Mn^^. 

In  this  section  of  the  report  we  have  shown  (i)  the  application  of  the 

new  formalism  enables  for  the  first  time  the  direct  calculation  of  Coulomb 

and  Exchange  crystal  field  parameters  and  (ii)  the  refined  theory  gives  a 

proper  interpretation  of  the  detailed  optical  spectra  unlike  the  simple  three- 

parameter  B,  C,  A  theory.  As  may  be  seen  from  our  published  work^,  the 

application  of  the  new  formalism  to  evaluate  the  Coulomb  and  exchange 

interaction  constants  would  require  solving  the  inverse-eigenvalue  problem 

2+ 

uisng  the  observed  optical  spectra.  The  observed  lines  for  MgF2:Co  from 
the  data  of  Sibley  ^  were  used  along  with  our  matrix  elements  and  the 
general  formula  derived  in  our  present  study  relating  the  diagonal  elements 
of  complementary  configurations. 

In  solving  the  inverse  eigenvalue  problem  and  obtain  V^'s  (i  =  1  to  10) 


it  is  necessary  to  have  a  sufficient  number  of  lines  from  optical  spectra. 

2+ 

For  MgF^.’Co  ,  eleven  well  defined  lines  were  available  and  ten  of  these  were 
used  in  the  inverse  eigenvalue  problem  and  the  remaining  one  was  used  to 
check  the  validity  of  the  solution.  To  our  knowledge  there  are  no  calculated 
values  of  Coulomb  and  exchange  interaction  constants  in  the  systems  we  have 
investigated.  Table  IV  gives  the  Coulomb  and  exchange  crystal  field 
parameters  and  a  comparison  with  the  corresponding  values  for  the  free  ion 
calculated  by  us  using  the  wavefunctions  given  by  Clementi  (see  Ref.  1). 

Also,  we  have  shown  in  Table  V  the  comparison  of  the  observed  and  calculated 

or 

energy  values  together  with  previous  tentative  assignments  and  the  revised 

assignments.  It  is  very  significant  that  our  analysis  has  provided  the 

correct  assignments  to  many  of  the  lines  as  shown  in  Table  V. 

The  application  of  the  inverse  eigenvalue  problem  to  the  data  on 
2+  5 

MgF2:Mn  with  our  matrix  elements  for  d  configuration  resulted  in  a  matrix 

2+ 

of  dimension  43x43  compared  to  20x20  for  Co  .  The  values  of  the  parameters 
for  Mn  case  are  also  given  in  Table  IV.  Again  we  have  compared  them  with 
corresponding  free  ion  values  calculated  by  using  dementi's  wavefunctions. 

As  was  pointed  out  earlier  there  are  no  previous  values  for  V.j's  from  other 
sources  for  comparison.  We  have  thus  shown  for  the  first  time  it  is 
possible  to  obtain  them  using  our  refined  formalism.  Also,  we  have  determined 
that  many  of  the  tentative  assignments  of  the  lines  belonging  to  this  system 
have  to  be  revised  (see  also  our  paper  Ref.  1). 


System 

V  -V 

V  -V 

V  -V 

''l  ''s 

^^3 

''6 

^^7 

''8 

''9 

'^10 

10  Dq 

MgF^rCo^'^ 

13034 

13728 

-2433 

2725 

8491 

3219 

1928 

1259 

6257 

-7250 

Free  Ion* 

16920 

16920 

0 

4379 

9724 

5932 

2189 

2189 

8460 

— 

MgF2:Mn^'^ 

11176 

11931 

647 

2469 

6026 

2725 

1584 

2164 

6093 

5710 

Free  Ion* 

15328 

15328 

0 

3958 

8806 

5378 

1979 

1979 

7664 

— 
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Table  V  Observed  and  Calculated  Energy  Values  in  Cm  for  MgE^^Co 


Previous 

Assignment* 

Present 

Assignment 

Observed  Energy 
Value 

Calculated 
Energy  Value 

\  ('f> 

7519 

7518.9 

% 

%  ('g) 

11905 

11904.9 

\  ('f) 

%9 

15385 

12205.1 

\  ('g) 

17452 

17452.1 

'Ti9  (*P) 

r 18519 

J  19305 

19305 

L 20921 

20833 

20832.9 

\  ih) 

21739 

21739 

%g'  '^ig  ('h) 

\  ('h) 

r23095 

23094.9 

'^ig  ('«> 

V23529 

23529 

S  ('h) 

'Eg  (^H) 

25000 

26968.9 

*  Ref.  26. 


6.  Conclusion 


(1)  In  this  program  we  have  examined  in  detail  the  photoluminescence  data 
on  InP:Fe  and  the  low  temperature  optical  absorption  results  on  GaAs:Fe  and 
derived  the  set  of  crystal  field  parameters  that  will  best  describe  the 
experimental  observations.  The  results  can  be  explained  as  intracenter 
transitions. 

(2)  For  the  d^  configurations  of  defect  centers  we  have  proposed  a  new 
refinement  for  the  theory  of  electron-electron  interactions  and  derived  the 
generalized  d-electron  matrices  for  d’^(n  =  2, 3, 4, 5).  Further,  we  have  derived 
the  significant  relation  connecting  the  complementary  configurations  d"^  and 

d  ^  ’  and  shown  that  in  the  generalized  case  the  diagonal  elements  of  the 

matrices  do  differ  by  varying  amounts.  We  have  applied  our  generalized  theory 
2+  3+ 

to  GaAsrCr  and  GaAs:Cr  cases  and  obtained  energies  in  good  agreement  with 
previous  work  by  Hemstreet  and  Dimmock  who  proposed  modifications  of  e  and  t^ 
orbitals  using  parameters  from  X  calculations. 

'I 

(3)  As  part  of  the  defect  center  problem,  we  have  examined  the  effect  of 

Jahn-Teller  distortion  and  presented  a  molecular  orbital  treatment  in  the 

2+  2+ 

framework  of  the  ligand-field  theory  for  the  cluster  Cr  -As^  in  GaAs:Cr  . 

The  charge  transfer  effects  have  been  derived  and  the  parameters  that  will 
interpret  the  EPR  results  of  Krebs  and  Stauss  have  been  obtained.  These 
results  are  significant  for  accurate  calculations  of  Coulomb  and  exchange 
interactions.  The  effects  of  charge  transfer  are  shown  to  be  very  important. 

(4)  We  have  also  verified  the  use  of  the  revised  formalism  for  d-electrons 
to  two  other  systems  (MgF2:Co  ,  MgF^tMn  )  and  derived  for  the  first  time  the 
complete  sets  of  crystal  field  parameters.  Also,  we  showed  the  advantage  of 
the  refined  treatment  to  interpret  correctly  and  in  detail  the  observed 
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Generalized  </*electro0  interaction  matrices-their  derivation  and  impact  on  existing  results 

R.  R.  StaAmia,  M.  H.  de  A.  Viccaro,  and  S.  Sundaram 
Department  cf  Physics,  Viiiversityc/ Illinois,  Chicago,  Illinois  60680 
(Rtoeived  2  Jtouuy  1980) 

The  ftstrictioii  of  the  “puie"  d  otbiuk  on  the  Tenebe-Suaeno  electron-electron  interaction  matrices  has  been 
removed  and  new  matrices  have  been  derived  appropriate  to  d-electron  transition-metal  ions  in  solids  and 
complexes.  We  have  adopted  the  group-theoretical  technique  developed  by  Racah  involving  the  coefficients  of 
fractional  parentage.  A  general  expression  useful  for  obtaining  the  matrix  elements  for  the  complementary  states  has 
also  been  given.  The  significance  of  the  present  study  has  been  emphasized  by  checking  and  comparing  our  results 
with  the  recently  published  results  on  QaAs£t^*  and  CaAs:Cr^*  and  by  presenting  the  examples  of  MgFjiCo^^ 
and  MgF,:Mn’'^  where  it  has  been  passible  to  deduce  the  values  of  the  d-electron  Coulomb  and  exchange 
parameters  in  conjunction  with  the  fine  experimental  optical  spectra.  Other  applications  such  as  the  removal  of 
accidental  degeneracies  have  been  indicated.  Important  points  for  further  improvements  have  been  discussed. 


I.  INTRODUCTION 

There  has  been  a  considerably  increasing  inter¬ 
est  in  the  study  of  the  transition- metal  ions  in  com¬ 
plexes  and  solids.  The  present  study  is  part  of  our 
systematic  investigations  on  the  defect  centers  of 
transition-metal  ions  in  III-V  semiconductors  of 
GaAs  and  InP.  Quite  recently  Hemstreet  and  Oim- 
mock*  have  calculated  many  electron-crystal-field 
term  states  ibesides  spin-polarized  spin  states)  for 
Ci*'^  and  Cr**  in  GaAs  using  an  improved  approach 
that  is  a  modification  of  the  standard  electron- 
electron  interaction  matrix  elements***  by  introduc¬ 
ing  parameters  accounting  for  the  changes  in  the 
single-particle  electronic  states  in  the  system. 

The  parameters  were  deduced  from  the  Xa  scat¬ 
tered-wave  cluster  calculations  which,  taken  alone, 
could  not  account  tor  the  electronic  structure  of 
Cr  in  GaAs,  In  the  process  of  our  investigations 
of  the  electronic  structures  of  Co**  and  Mn**  in 
solids,*  and  in  our  attempts  to  compare  our  results 
for  Cr  in  GaAs  with  those  of  Ref.  1,  we  have  en¬ 
countered  the  difficulties  associated  with  the  un¬ 
availability  of  the  explicit  generalized  form  of  the 
Tanabe-Sugano  matrix  elements  for  the  single¬ 
particle  electronic  wave  functions  which  are  not  of 
the  “pure”  d  character. 

While  Hemstreet  and  Dimmock  have  suggested 
the  modified  approach  of  removing  the  restriction 
of  pure  d  nature  of  the  electronic  wave  functions  by 
introducing  three  parameters  (,R„,  and  R,,), 
the  advantages  of  a  generalized  approach  have 
prompted  us  to  derive  the  explicit  relations  and 
provide  a  general  formalism.  We  therefore  pre¬ 
sent  in  this  paper  the  explicit  and  general  forms 
of  the  Tanabe-Sugano  matrix  elements  with  the  aim 
of  making  them  easily  available  to  the  future  re¬ 
search  workers.  We  have  also  used  them  here  to 
check  and  compare  our  results  with  the  calcula- 
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tions  of  Hemstreet  and  Dimmock*  for  Ci**  and  Cr** 
in  GaAs.  As  will  be  shown  in  Sec.  HI,  there  is  a 
substantially  good  agreement  between  their  results 
(Ref.  1)  and  those  obtained  from  our  present  gen¬ 
eralized  approach.  The  general  matrix  elements 
are  also  useful  for  making  refined  calculations  of 
the  transition-metal  ions  in  solids.  Recently  we 
have  utilized  them  for  designating  correctly  the  en¬ 
ergy  levels  and  for  deducing  the  Coulomb  and  ex¬ 
change  parameters*  as  well  as  the  refined  values 
of  the  cubic-field  splitting  parameter  associated 
with  Co?*  and  Mr?*  ions  in  M^,.  Moreover,  the 
general  treatment  is  helpful  for  removing  the  ac¬ 
cidental  degeneracy*  observed  in  certain  areas. 

In  Sec.  n  we  give  the  derived  general  results  with 
only  a  brief  outline  of  the  method  since  the  method 
is  already  well  known  .***■*  Section  HI  presents  ex¬ 
amples  which  emphasize  the  effects  of  the  general 
expressions  on  the  existing  results  in  the  literature 
to  justify  the  importance  of  our  derivations.  Sec¬ 
tion  rv  deals  with  the  discussion  and  conclusion. 

II.  BRIEF  OUTLINE  OF  THE  THEORY  AND 
DERIVATION  OF  MATRIX  ELEMENTS 

The  lack  of  knowledge  of  the  correct  wave  func¬ 
tions  for  the  transition- metal  ions  in  a  complex  or 
a  solid  has  been  a  stumbling  block  to  real  progress 
in  the  theoretical  understanding  of  the  experiment¬ 
al  observations.  Finkelstein  and  Van  Vleck*  were 
the  first  ones  who  treated  the  complex  ions  in  crys¬ 
tals  by  perturbing  the  free-ion  energy  levels  due 
to  the  crystal  environment.  Later  studies  by  many 
authors**’ •*•*■**  have  helped  further  our  knowledge 
in  this  field 

Racah*  developed  a  theory  using  coefficient  of 
fractional  parentage  to  explain  the  spectra  of  free 
ions  and  introduced  three  basic  electron  Coulomb 
repulsion  parameters  A,  B,  and  C  for  the  pure  d 
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electrons.  This  theory  was  extended  to  the  case  of 
transition-metal  ions  in  solids  and  complexes  by 
incorporating  the  effect  of  the  crystal  field  in  terms 
of  the  splitting  parameter  A  with  the  assumption 
that  the  d  electrons  retain  their  pure  d  character. 
This  theory  has  been  widely  adopted  for  interpret¬ 
ing  the  optical  spectra  of  transition  metal  ions  in 
solids  and  complexes.  In  fact,  the  theory  considers 
only  three  parameters,  B,  C,  and  A,  as  the  pa¬ 
rameter  A  plays  no  role  in  accounting  for  the  ener¬ 
gy  differences.  While  this  theory  brings  forth  use¬ 
ful  simplifications,  it  suffers  from  the  serious 
drawbacks  of  neglecting  the  important  solid-state 
effects,  such  as  the  modifications  of  the  electron 
wave  functions  from  the  pure  d  character.  The 
faults  of  the  theory  are  even  more  noticeable  when 
one  recognizes  that  it  really  does  not  work  well  in 
many  cases.*'*^**  For  example,  a  poor  agreement 
has  been  found  by  Stevens”  in  the  case  of  Co  com¬ 
plexes.  Also,  the  experimental  observations'*  of 
the  optical  absorption  of  Mii**  complex  reveals  a 
fine  structure,  consisting  of  two  peaks  separated 
1^  300  cm*',  which  cannot  be  accounted  for  by  con¬ 
sidering  spin-orbit  interaction  or  departures  from 
cubic  symmetry.  This  splitting,  however,  can  be 
easily  e^lained  on  the  basis  of  the  present  im¬ 
proved  theory  in  which  the  outer  electrons  are  not 
restricted  to  possess  the  pure  d  character.  In 
that  circumstwce,  the  accidental  degeneracy  in 
*Ai,  and  predicted  by  the  B,  C,  A  theory  is  in¬ 
nately  removed.  The  mechanism  suggested  by 
Koide  and  Pryce”  in  terras  of  the  mixing  of  the  odn 
lattice  vibrations  to  explain  the  *Ai,-*E,  splitting 
in  complexes  is  really  an  additional  effect  be¬ 
sides  the  one  mentioned  above. 

Thus  it  is  imperative  to  improve  the  theory  by 
removing  the  restriction  of  the  pure  d  nature  from 
the  electronic  wave  functions.  To  this  end,  we 
conceive  of  the  one- electron  orbitals  as  sets  of 
orbitals  which  form  bases  for  irreducible  repre¬ 
sentations  of  the  symmetry  group  of  the  crystal 
potential.  We  shall  follow  the  well-known  strong- 


TABLE  L  List  of  Independent  Coulomb  and  exchange 
Integrals. 


Alternate 

notettons 

Integral 

a 

<4fl«Vr„Uf) 

Vi 

b 

Vi 

e 

<94le’/r„|«4) 

Vl 

d 

<«4leVr„|€4; 

Vi 

e 

V% 

f 

Vt 

g 

(9fl|e*/r,ili7il> 

Vt 

k 

<««|s*/robi9) 

V% 

i 

) 

(44ls*/r,iM; 

field  coupling  scheme*-*  since  many  physical  cases 
of  interest  fall  in  this  category.  In  this  scheme, 
Tanabe  and  Sugano*  have  derived  the  <f- electron  in¬ 
teraction  matrices  in  the  simplified  picture  of  pure 
d  electrons  by  utilizing  Racah’s  irreducible  tensor 
operator  formalism  appropriate  to  octahedral  sym¬ 
metry.  We  have  adopted  the  same  formalism  and 
have  derived  the  improved  matrix  elements  in 
terms  of  the  ten  independent  Coulomb  and  exchange 
parameters  (Ref.  3)  V,,  V,, . . . ,  which  have 

been  tabulated  in  Table  I  in  standard  Dirac  notation 
for  easy  reference.  The  K/s  are  the  two-electron 
matrix  elements  of  the  Interaction  operator 
between  two  electrons  separated  by  the  distance 
ri2.  The  functions  B,  e,  i),  and  have  been  used 
to  denote  one-electron  orbitals  in  the  cubic-field 
representation.  Evidently  these  orbitals  are  no 
longer  of  the  pure  d  character. 

The  method  of  derivation  of  the  electron- electron 
interaction  matrices  is  well  known*  and  straight¬ 
forward  but  very  tedious  and  time  consuming.  It  is 
not  planned  to  give  details  here  because  it  has  been 
described  in  many  places.  The  electron-electron 
matrix  elements  for  d'  electrons  (with  «  =  2,  3,  4, 
and  5)  derived  in  the  present  study  are  listed  in 
Table  II-V  for  octahedral  symmetry.  It  should  be 

TABLE  II.  Tabulation  of  the  electron-electron  matrix 
elements  for  <st*  electrons  in  terms  of  the  Coulomb  and 
exchange  integrala  F/a  (a,  ft . etc.).  Asterisk  de¬ 

notes  the  change  in  sign  with  respect  to  Tanabe-Sugano 
matrix  elements  because  of  the  phase  conventions. 


‘A,(d*) 

t|('A,) 

e*(U,) 

e+/ 

•J6g+/2h 

a  *2} 

‘£(d*) 

e*(’B) 

f|(‘£) 

e*(‘£) 

e-f 

Zh* 

Zh* 

a-j 

‘TjW*) 

t|(‘r,) 

tjl’T, )«(*£) 

til'Tj) 

b*} 

-2<* 

tj(*T,)e(’B) 

-2f* 

d  +g+i/3il 

-c/V5 

’T,(d*) 

f|(*r,) 

fi« 

»|(’T,) 

b-i 

-z-rSi* 

tie 

-2/3i(* 

d  — 

’Aj(c*)=e-3/ 

’T,(sf2)  a/>/Dc-V5* 

‘r,(e<,)  =<f  +g+i/3i:-  (1/V3)h 
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noted  that  we  have  followed  the  phase  conventions 
of  Griffith*  in  our  derivations.  As  expected,  our 
matrix  elements  reduce  to  the  Tanabe-Sugano  ma¬ 
trix  elements*  in  case  one  takes  the  limiting  condi¬ 
tion  that  the  0,  e,  t},  and  £  orbitals  are  pure  d- 


type  orbitals.  Those  matrix  elements  in  Tables 
II -V  which  differ  in  sign  on  reduction  to  pure  d- 
orbital  case  from  the  corresponding  matrix  ele¬ 
ments  obtained  by  Tanabe  and  Sugano*  (or  by  Griffith*) 
have  been  marked  (in  the  tables)  by  an  asterisk. 


TABLE  in.  Electron-electron  matrix  elements  for  d*  electrons.  Asterisk  denotes  the 
change  in  sign  with  respect  to  Griffith’s  matrix  elements. 


/J(*£)e'‘(‘A,)  f|(‘A,)e‘(*£) 

f|{'£)e‘(*£) 

f“('A,)e’(*£) 

36  -2^/6t 

-n/«* 

0 

/i(U,)e‘(*£) 

-2/6i  a  +  (2//3)c+2d 

-g-(l/>/5)6+2; 

(4/v^c+(2//3)6*  v/3^  +  6 

-V^’  4c//3  +  (2/v'3)** 

a  +(2//3)c+2(f 

2h* 

0  'fZg  +  h 

2h* 

3e-5/ 

tteh^Ai) 

f2C*(‘£) 

f|  j+26- 

■2j  -Vli*  vTi* 

0 

-2h 

f|(’T,)e 

6  +  (4//3)c+2d  ->/3* 

3i* 

b*2d-g 
-<2//3)A  +;• 

-^i 

-f* 

(2/>/l)c+2d  +  c 
-3/  +i'  +  (l//3)6 

+2h* 

fj«*(‘£) 

e+2d  +  (2/V3)c 
-^-(6 //§)-/ 

%  <1 

<i(’T,)e  <|(‘r2)e 

tieH'E) 

t|  a +26  — 3i  — 5t* 

2g  +  26//3 

2h//3 

fi<*r,)e 

b-}+2d+g  /3h* 

+2A//3 

-37 

-3: 

f|(‘rj)e 

b*4c/'f3*2d 

~g*} 

-i* 

+7* 

e  +  2d  +  2c//3 
-g*f-h/^ 

4c/V3+26//3 

2c/V3+2d  +  e-/ 

— g  —  h  //3 

*T, 

f|<"r,)e 

tieH^At) 

fltV,)* 

6+4c/Vl+2d 

-2£-> 

2/3i 

2/» 

2cA^  +  2d  +  c 
-3/  -  2^ 

-26//3 

<i(‘£)e:  *A,-a-2cA^+2d 

-g-/3k-j 

ff('£)e:  *A,  a+2v^  +  2d- 
tl:  ‘a  2 -36 -37 

■  g+h//3  —  } 

%  =  b*2d-2g-4hAf3-j 

V-’ff) 


If 

li(*T,)e 


a-^So  — 3/  —  \2i 


-VDf 


TJ"  -TT" 


a  +  26  +  v'Sc  2c  +  6 

*3d-2g 
-  V3*-2i 


3  +  26+ %'Jc  -\^«*  ->/J<* 

*Zd-2g 


4  .  4  2. 


+  e+/-2/f 

-h"-’ 


b^‘^^c-¥Ad  —26 


+  e-/-2^ 

2  u 

"VJ*"-' 


6  +  -^c  +  4rf  v/^» 

i/5 

+  e_3/-2^ 

1 

•^c  +  3(^  +  3c 


ri<»r,)e  /i(*T,)c  /|f7',)e*<l4j)  )«*(>£)  t|(‘rj)c*(‘i4i)  (jC> 


4  a*5b.j 

V6j  -5V2« 

u 

-75" 

u 

0  +  26+ V3c  -2c -6* 

VTi 

-/3« 

— /5« 

-/26* 

*3d—'/Sh~2j 

o  +  26  +  V5c 

-3«* 

+5«* 

-5i* 

/5y+^6 

+  3rf+  \^6 

b-¥-^c*id  +  e 
-2f+2g*-^h 

-2/56 

0 

-3v'5i* 

-j 

4  2 

-7?" -75* 

4(1  r, )«*<«£) 

6+-^c+4rf+e 

2 

+/5<* 

.f-2g.^k 

*} 


4^T,)«*(‘i4,) 

*■*’75'' 

+/5<* 

+  e  +/-2^ 

5 

<L-* 

-7*c  +Srf+  3« 
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TABLE  IV.  (Continued) 


‘Aj  t\  tieE)e  /|(‘E)c'2(‘£) 


4 

4^E)e 

/!|(‘Ai)p2(‘Ai) 


l|('£)c*('£) 


2fl  +  46  — 4\7>j 

36  +  \'7c  +  3rf 
+  2x'36 


-4v75i 


4 

a+-;sc  +  4d 

\  O 

*e  +  f-2g 

-^/,  +  27 


n  +-?5C  +  4d  *t’ 
'  ^  2 
-/-2a-- -4/1 


■7?" 


0 

0 


v'(rA'+ 1"?* 


+  2%'^/i* 


6(’-10/' 


*£ 


<?(2£)e  I|(‘£)^2(‘A,)  4(‘A,)e2(‘£)  <§(»£)<-*(»£) 


4 

2a  +  46  -  3j  2  \^i 

/*  » * 

4fE)e 

36  +  3rf 

— 4\/7» 

4(‘£)e*(‘Ai) 

4 

a  +‘^c  +  4rf 

4\/J  2v'5,, 

+  e+/-2/f 

4(‘A,)«*(‘£) 

2  , 

-75"-' 

a  +“^c  +4rf+  <» 
2 

-f-2g~-ysh 
+  2> 

4(‘£)e»(‘£:) 

0 


4^■2■  2v'2, 


0 


4 

a  +  "TBpr  +  4</ 

^  2 
+  e-/'_2/? 


4(®r,)e 

a  +  2h+  /Sc  +  3d  2c  +  *  >/^»* 

+  /5h  -  2} 

/Sh 

4<'T’j)c 

rt  +  2ft  +  /Sc  — 5f  * 

+  3d—/Sh  ^ 

4(‘li)e* 

6+_c  +  4rf 

-/Si" 

+  e-f-2g 

2  ^ 

-75 

r  4 

-^*3d-<r3e-3f—j^h 

’Tj 

4(*r,)e 

t\^T{)e 

4(’T,)c*(’a,) 

tl(^T,)e^(‘£) 

h<^ 

a  +  2ft  +  V7c 

« 

1 

C3 

1 

/Si 

.* 

f 

-,Sh* 

4frj)c 

/f(»7’,)e*(»A,) 

*Sd-2g-^ 

-2} 

a  +  2ft  +  /Sc 
+  3d -2^-  -JSh 

-/Si* 

4 

6  +  ^5?C  +  4</ 
v3 

2v/5  * 

-TT" 

_  2/2ft 

-2^'5i* 

4(>T, )?*(»£) 

+  e  —  3f-i 

6+^c  +  4d 

<,e» 

+  e-f-2g 

2  t  • 

-75"-' 

-^c  +  3d  +  3c  —  5f-2g—-^li 

te 


fE)e  tl^£)e‘fA,) 


1 


I 


4(*^)e  36  +  v^c  +  3rf  —*^^6  0 

+  Jf +  ■^-3/ 

^(*£)c  36+\/^+3rf  -VUi* 

<|(‘B)«*(l4j)  a  +  -^c  +  4rf 

+  e-3f-2g  --^h-i 

<J(*£)e  36+  y/^c  — 4^'^i 

+  3rf-2jf 

<i+-^c  +  4<# 

'  ■*  2 
+  e  _  3/—  2g'  —  -^h  +  2; 


(ff£V  3h  +  ^f^c  +  3d  2x'3»* 
-2%'J6 

<|(‘£)e*<‘£)  a+-^c  +  4d 


+  p-f-2g 


E)c  36  +  v^c  +  3(i 

-2g-Ah/ 

j|(‘i4j)e  ^E  *  3h  ■>■  yf^c  ^  3d 
—  3g—  \^h  —  3i 
«|(>r,)e*(*.4j)  %»b  +  4c/y,'^*4d 
+  e-3f—4g 
-Ah/yf^-j 


Although  we  have  followed  consistently  Griffith’s 
phase  conventions,  we  have  not  been  able  to  trace 
out  why  we  disagree  in  sign  in  some  of  our  ma¬ 
trix  elements  from  those  listed  by  Griffith.^ 

As  for  the  matrix  elements  of  the  complementary 
states  <!“*■•  («  =  2,  3, 4,  5),  the  nondiagonal  matrix 
_ I 


elements  remain  the  same  as  in  d"  matrices  but 
the  diagonal  matrix  elements  of  d*®""  states  differ 
from  the  corresponding  ones  of  d"  by  a  function 
vdiich  is  a  linear  combination  of  V,'b.  The 
straightforward  but  tedious  group-theoretical 
method  yields  the  general  expression, 


< /S'  "(s,r,  k*  • '  )srA/s  Mr  |  (eVrJItS'”  (s^r,  '(^r*)srA/s  Af  p) 

-  <  k '  )SrMs  |  (e»/r„ )  |  / J'(S,ri)e'{S4r„  )SrM^Mr} 

=  (3  -  m)Vt  +  4(3  -  w)V,  +  (24  -  4w  -  6/)(F,/'/7  +  +  (6  -  3l)V, 

+  (5/-10)Ve  +  (2w  +  3/-12)F,  +  (2»i  +  3/-  12)7*/ /J- 2(3 ,  (1) 

where  the  integers  m  and  I  are  connected  to  n  by  the  relation 

m  +  /  =  «  or  10  -  «.  (2) 

InEq.  (1)  <2  represents  the  orbitals  and  £  and  e  represents  the  orbitals  0  and  e.  S,  and  S  are  the  tot¬ 
al  spin  quantum  numbers  and  r/s  and  r  are  the  irreducible  representations  of  the  cubic  group.  The  dif¬ 
ference  between  the  diagonal  element  of  d*  and  d”'"  as  given  by  Eq.  (1)  reduces,  as  expected,  for  the  pure 
d-electron  case  to 


<  /2*'"(Sir,)^''(5ir2)SrAf,Mr  |(eVru)|tJ-"'(S,r2)e«''(S»ra)SrA/,A/r> 

- </?(S,r>'(5arj)SrAfsAfr  |  (e^/ru)!  t;'(Stri)e'(S,rt)SrMgAfr)  =  (45A  -  70B  +  35C) -  (9A -  14B  +  7C)(w  +  1) , 

(3) 

I - 


where  A,  B,  and  C  are  the  standard  Racah  param¬ 
eters.  One  notes  that  in  the  case  of  pure  d  orbit¬ 
als,  the  expression  (3)  depends  only  on  the  sum 
(m  +  /)  which  according  to  Eq.  (2)  is  n  or  10  -  n  and 
hence,  for  a  given  n  (with  whatever  appropriate 
values  of  the  set  w,  1)  the  difference  between  the 
diagonal  elements  in  Eq.  (3)  is  constant.  The  above 
formula  is  particularly  significant  since  it  reveals 
that  the  diagonal  elements  of  the  d*®*”  configuration 
change  by  different  amounts  from  the  correspond¬ 
ing  ones  for  the  d”  configuration.  As  an  example, 
for  the  configuration  d®  the  diagonal  element 


<tJ(*£)e®(*£)»E|(eV»'u)U®(®£)c®(‘£)»£) 

and  the  corresponding  one  in  the  configuration  d 
differ  by  the  amount 

[-  12(7,/'/7  +  K.)  -  67*  + 107*  +  67,  +2/J7.] , 

whereas  the  diagonal  element 

<t»(*A,V»(*£)*£|(cVru)U»(‘A,)e»(»£)*£> 

in  the  configuration  d’  differs  from  the  corres¬ 
ponding  one  in  the  d^  configuration  by 
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TABLE  V.  ElectroD-«lectroa  matrix  elements  for  electrons.  Asterisk  denotes  the  change  in  sign  with  respect 
to  Griffith’s  matrix  elements. 


2b  +  43  + -^c +4<f  -  2g- 3> 

— 1/&'  0 

3b+2/3c  +  6d  -4k* 

*e-/-3g 

-/3h 

36+2/3c+6d 

0 

*e-3f  +  3g*-^h-3j 

f|(‘E)«’ 

*A,  t\0E)e 

f|(*£)e*(‘E)  t|(‘E)c* 

10 

a  +  +  6rf 

+3e-3f-3g--^pi-j 

<!('£)* 

2a+4b*4d 

g->f3k* 

-2£-^* 

-3} 

f|(*£)«*(‘£) 

36  +  2/3c  +  6d 

-VW* 

*e-f-3g 

-yfSh 

fi('£)«’ 

2 

a  +  *  ®d 

+3c-5/-3g 

I 


I 


’£ 

<}('£)« 

f|(’£)B’(U,) 

f|(*£)e*(’A,) 

f|(*£)e*(‘£) 

fi('£)c’ 

fi(‘A,)e’ 

2«*4»-^ 

7r*7T’‘ 

2/f< 

6i 

2/6i 

fl(‘£)* 

*4d-2g-^ 

2a  *4b  +  -^c 
+4d-2g 

-m 

3< 

0 

fK*£)**(U,) 

f^*£)BVAt) 

-kH.3) 

33+2>f3c+6d  0 
+  ,+/-3g 

33+2/3c 

+6d+e-3/-g 

0 

ZVIk 

1  1 

• 

-Z/3i 

6i 

fi(*£)#*('£) 

If 

-7f 

3b+Z^/lc+6d 
+  t-/-3g 
-/3k 

0 

2/li 

fi('£)** 

a  +2/3c  +G(i 

4  2  . 

-TT-TT* 

+  3e-5/-3g 


a  +  2^/3c*Gd 
*3e-5f-3g 


‘T, 


36  +  2i/3c+6<i  -2h 

*e-3/-4g 


3b+2/3c+6d 
+  e-J-3g 
-/3h~3} 


*Ti  t\(^Tt)e 


#4/3-7 


tl{%)eH^Ai) 


2 

a  +  5b  ^  ~7^c 
V  3 

+  4({-3g-\f3h 
-3j 


—  l/W* 

g--^>r 

a  *  26  ■+  2v  3f 

—  \  6f 

-■ig-  -Y  ^ 

*  4  -—c  ^  Gd 

V  3 

3e  —  5y  4^ 

a+56  +  2V3c 
+  4<i  —  3^— V3// 
-3> 

/2(^7’,)e^l42) 


i/2/* 


a  +  2b  +  2/3c 
+  6d*e-3f 


-2} 


g  +\'3)i* 


-V2,- 


6  +  ^  c  3d 
V  3 

+  3e  -  5/  -  ig 
4  ,  . 

-  ^A-; 


[-3K,  -  12V,  -  6(V3/>^  +  V4)  +  3V,  -  5V, 

+  3V2+V7Vg  +  6Vio]. 

Also,  in  the  configuration  the  diagonal  element 

{tl{^T,)e^WT,HeVr„)\tlWeH^EVT,) 

differs  from  the  corresponding  diagonal  element  in 
the  (f'  configuration  by 

[-  V,  -  4  V,  -  10(  V,/V7  +  VJ  -  3  V,  +  5  K, 

+  5V2  +  5V,//5  +2VjJ  . 

It  is  clear  from  the  above  examples  that  the  vari¬ 
ous  diagonal  elements  are  changed  by  the  different 
amounts  in  going  from  the  d’  to  configuration. 
On  the  other  hand,  if  pure  d  oribitals  are  involved, 
all  the  diagonal  elements  of  the  d®  configuration 
differ  from  the  corresponding  diagonal  elements  of 
the  d'  configuration  by  the  constant  amount  (-18A 
+  28B-  14C)  which  agrees  with  the  result  given  by 
Griffith.*  Thus,  in  general  the  present  results 
from  our  improved  treatment  differ  from  the  cor¬ 
responding  results  in  the  simplified  theory  of  the 
pure  d  orbitals  where  a  constant  difference  has 
been  well  known.  The  significance  of  the  new  di¬ 
agonal  elements  will  be  made  clear  in  Sec.  HI. 

It  must  be  recalled  that  in  cubic-crystal  fields 
the  diagonal  matrix  elements  for  the  states 


are  further  admixed,  as  usual,  by  the  crystal  field 
term  (0.4^-  0.6p)A,  where  A  =  iODq  is  the  cubic- 
field  splitting  parameter.  The  splitting  parameter 
A  is  positive  for  the  octahedral  symmetry  and  in 
that  case  the  level  lies  below  the  e  level.  On  the 
other  hand,  in  the  case  of  tetrahedral  symmetry  it 
has  been  explained  by  Griffith*  that  the  parameter 
A  changes  sign  thereby  inverting  and  e  levels. 
Also  the  diagonal  matrix  elements  of  the  d"  con¬ 
figuration  are  related  to  those  of  d*®'"  configura¬ 
tion  in  a  given  symmetry  (octahedral  or  tetra¬ 
hedral).  This  is  a  very  important  result  concern¬ 
ing  our  improved  theory  since  the  diagonal  ele¬ 
ments  no  longer  change  by  the  same  amount  as  ex¬ 
pressed  mathematically  by  Eq.  (1).  Thus,  in  sev¬ 
eral  cases,  based  solely  on  this  type  of  change  in 
the  diagonal  elements,  the  energy  levels  in  the  new 
treatment  are  expected  to  be  different  and  certain 
levels  which  are  degenerate  in  the  old  treatment 
are  now  expected  to  have  split  components  in  the 
improved  new  treatment. 

111.  IMPACT  OF  IMPROVED  MATRICES  ON 
EXISTING  RESULTS 

At  first  glance,  the  full  import  of  the  matrices 
derived  in  Sec.  II  may  not  be  obvious  owing  to 
the  fact  that  one  now  has  eleven  parameters 
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1^2 . and  a)  replacing  four  parameters 

{A,  B,  C,  and  A)  of  the  old  theory.  The  large 
number  of  parameters  makes  the  theory  apparently 
less  attractive,  but  it  is  important  to  note  that  it 
includes  the  effect  of  the  surroundings  in  the  cor¬ 
rect  way.  The  new  theory,  however,  is  very  sig¬ 
nificant  since  it  has  been  able  to  (i)  test  the  calcu¬ 
lations*  made  on  GaAs ;  Cr**  and  GaAs :  Cr**  and 
make  comparisons,  (ii)  give  correct  group-theo- 
retical  assignments  to  the  observed  energy  levels 
of  the  impurity  in  solids  and  predict  the  values  of 
the  Coulomb  and  exchange  interaction  constants 
which  are  otherwise  not  possible  to  know  from  any 
other  source,  and  (iii)  remove  the  accidental  de¬ 
generacy  inherent  in  the  J3,  C,  A  theory.  It  is  con¬ 
ceivable  that,  in  due  course,  from  the  applications 
of  the  present  matrices  to  a  variety  of  physical 
problems,  the  full  significance  of  the  theory  and 
the  generality  of  the  present  superior  formalism 
will  become  even  more  obvious.  Only  the  first  two 
points  will  be  discussed  in  the  following  section 
since  the  third  point  has  already  been  discussed 
briefly  in  Secs.  I  and  II. 

A.  GaAs:Cr**'  and  GaAs:Cr**  energies 

The  impurities  Ci^*3d*  and  Cr**3r/*  occupy  the 
tetrahedral  site  substituting  Ga  in  GaAs.  As  men¬ 
tioned  in  Sec.  I,  Hemstreet  and  Dimmock*  have  re¬ 
cently  performed  calculations  for  these  systems  by 
modifying  the  free- ion  one-electron  orbitals  of  fj 
and  e  symmetryinsuchawaythattheir  normaliza¬ 
tion  constants  were  altered  by  means  of  the  pa¬ 
rameters  R„,  R„,  and  R,,  which  were  deduced 
from  the  xa  calculations.  We  have  used  the  matrix 
elements  listed  in  Tables  in  and  1/  Jpr  d"'  and  d* 
electrons  to  verify  the  earlier  calculations*  for 
GaAs :  Cr**  and  GaAs :  Cr**.  Repetition  of  the  cal¬ 
culations  yields  general  agreement  of  our  results 
with  those  of  Hemstreet  and  Dimmock.*  However, 
there  are  some  minor  differences  [see  Fig.  (I)]: 

(i)  their  two  *Ti  levels  of  energies  0.74  and  0.84 
eV  give  energies  0.80  and  1.02  eV  in  our  calcula¬ 
tions,  (ii)  ^Ti  in  Ref.  1  should  have  been  designated 
’Tj,  and  (iii)  *£  with  energy  1.45  eV  and  ’Tj  with 
energy  1.33  eV  in  Ref.  1  should  have  been  origin¬ 
ated  from  the  configuration  instead  of  from  the 
configuration  t*.  The  above  slight  differences  may 
be  due  to  some  errors  in  their  d*  matrix  elements. 

We  have  also  calculated  the  energy  level  for 
GaAs :  Cr®*  and  compared  our  values  with  those  ob¬ 
tained  earlier.*  Figure  2  depicts  our  calculations 
where  comparison  has  also  been  made  with  the  re¬ 
sults  of  Ref.  1  (Keeping  the  same  values  of  the  pa¬ 
rameters  involved  as  in  Ref.  1).  Our  results  [Fig. 
2(b)]  are  in  good  agreement  with  earlier  results 
(Fig.  2(a)]  except  that  their  at  0.7  eV  and 
at  0.95  eV  should  now  both  be  *Ti. 


B.  Refined  interpretation  of  the  experimental  spectra 
and  deduction  of  exchange  and  Coulomb-interaction 
constants 

A  short  report*  on  the  application  of  the  present 
refined  treatment  to  the  analysis  of  the  recently 
observed*®  fine  optical  spectra  of  MgFj :  Co“*  d''  and 
R^Fj  :  Mn**rf®  has  already  been  presented.  Here  we 
only  briefly  state  the  concluding  remarks  in  order 
to  further  justify  the  importance  of  the  refined 
matrices  given  in  Tables  II -V. 

Sibley  and  co-workers*®  and  others  have  per¬ 
formed  excellent  optical  experiments  which  yield 
many  sharp  and  abundant  lines  of  the  transition- 
metal  ions  in  solids.  They  have  observed  not  only 
the  spin  quartet  states  but  also  many  hitherto  un¬ 
observed  spin  doublets  by  lifting  the  spin  forbid¬ 
denness  of  the  excitations  through  exchange  inter¬ 
actions  with  the  color  centers.  These  abundant  and 
fine  lines  offer  an  opportunity  for  their  interpreta¬ 
tion  by  means  of  the  refined  treatment  presented 
here  and  also  at  the  same  time  obtain  the  values  of 
the  parameters  V,’s  along  with  the  more  accurate 
value  of  A  by  solving  the  inverse  eigenvalue  prob¬ 
lem.  We  find  that  in  the  case  of  MgFj  :C(f*  [using 
Table  III  and  Eq.  (1)]  some  of  the  previous  tenta¬ 
tive  assignments*®  of  the  lines  based  on  the  three- 
parameter  theory  are  now  altered.  The  line  pre¬ 
viously  designated  as  *£,(*G)  is  found  to  be  really 
a  mixture  of  *£,(®G)  and  *3\j(®G);  *rjj(®G)  is  really 
“Tjj  (®G).  The  deduced  values  of  the  interaction  param¬ 
eters  and  the  splitting  parameter  AineV  are:  Fj 
-^2  =  1.616,  =  1.702,  1^1-^5=  0.302, 

=  0.338,  ^8=  1.053,  77=0.399,  78  =  0.239,  7, 

=  0.156,  7j8=0.773,  andA  =  -0.899.  Compared  to 
ttiese  values,  our  calculated  free- ion  values  (using 
dementi’s  wave  functions*®)  are:  7i  -  72  =  2.097, 
7,-78=2.097,  7i-  75  =  0,  73  =  0.543  ,  78=1.205, 

7,  =  0.735,  7b  =  0.271,  79  =  0.271,  and  7,o  =  1,049  eV. 

As  for  Ml*'*  3d®  in  MgF2,  for  the  second,  sixth, 
eighth,  and  ninth  lines,  the  assignments  given  in 
Ref.  15  as*ri,(*G),  *r„(*P),  ■'r„(*F).  and*r2,(*F) 
turn  out  to  be  ®r2,(®Z)).  Mi,(“G),  and 

*F),  respectively.  The  deduced  values  of  ex¬ 
change  and  Coulomb  parameters  and  the  refined 
value  of  the  cubic- field  splitting  parameter  in  this 
case  are:  7,  -  72  =  1.385,  7,  -  7^  =1.479,  7,  -  75 
=  0.080,  73  =  0.306,  78  =  0.747  ,  77  =  0.338,  79 
=  0.196,  79  =  0.268,  7,o  =  0.755,  and  A  =  0.708  eV, 
which  compare  with  our  calculated  free- ion  values 
(using  dementi’s  wave  functions*®)  of  7,  -  73 
=  1,900,  7i-  7,  =  1.900,  7i-  75  =  0,  73=0.491,  79 
=  1.092  ,  77  =  0.667,  78  =  0.245  ,  79  =  0.245,  and  7,o 
=  0.950  eV. 

Thus,  the  refined  matrix  elements  are  useful  to 
analyze  the  spectra  of  ions  in  solids  particularly 
when  abundant  and  sharp  lines  are  available  from 
experiments.  The  analysis  yields  the  Coulomb  and 


23 


c;  K  N  K  R  AL  I  /  K  I)  </  ■  K  I,  K  C  T  R  O  ^  I  T  F.  R  A  C  T  I  0  .%  M  A  T  R  I  C  K  S  .  .  . 


exchange  integrals  along  with  the  refined  values  of 
A.  Moreover,  the  symmetry  assignments  of  lines 
are  deduced  and  these  are  expected  to  be  correct 
since  the  treatment  is  based  entirely  on  group  the¬ 
ory.  Also,  this  analysis  will  stimulate  further  ex- 


Tin 

periments  to  observe  more  and  sharper  experi¬ 
mental  lines  and  also  further  encourage  the  devel¬ 
opment  of  the  first-principles  theory  to  explain  the 
deduced  interaction  and  crystal-field  splitting  pa¬ 
rameters. 


(a)  (b) 


FIG.  1.  Shows  the  energy  levels  of  GaA8;Cr^*  (a)  calculations  by  Hemstreet  and  Dimmock  (Ref.  1)  and  (b)  present  cal¬ 
culations  employing  d*  matrices  of  octahedral  symmetry  from  Table  IV.  The  values  of  the  parameters  A,  B,  C,  A, 

Hit,  and  R^,  in  both  cases  are  kept  the  same  as  in  Ref.  1. 
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IV.  DISCUSSION  AND  CONCLUSION 

The  need  to  improve  and  generalize  the  Tanabe- 
Sugano^  interaction  matrices  has  been  justified  in 
Secs.  I  and  m  in  order  to  eliminate  the  defects 
present  in  the  previous  treatment  which  is  based 
on  the  consideration  that  the  <f-electron  transition- 
metal  ions  in  solids  contain  pure  d  orbitals.  The 
improved  treatment  in  Sec.  11  incorporates  the 
solid-state  effects  adopting  group-theoretical  tech¬ 
niques  and  is  useful  for  analyzing  the  experimental 


results  more  accurately  than  before,  particularly 
when  abundant  and  fine  lines  are  available  through 
experiments.  Such  analysis  in  the  present  investi¬ 
gation  has  provided  not  only  the  definite  group- 
theoretical  assignments  of  the  observed  lines  but 
also  the  (hitherto  unknown)  values  of  the  Coulomb 
and  exchange  interaction  constants  and  refined  val¬ 
ue  of  A.  Also,  the  present  treatment  examines  and 
verifies  the  existing  calculations'  (see  Sec.  Ill)  on 
GaAs :  Cr**  and  GaAs :  Cr'*.  Furthermore,  the  im¬ 
proved  matrix  elements  (Tables  II-V)  can  be  used 
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FIG.  2.  Energy  levels  of  GaAs  iCr**;  <a)  Hemstreet  and  Dlmmock’s  calculations  (Ref.  1)  and  (b)  present  calculations 
using  d*  matrices  of  octahedral  symmetry  from  Table  IIL  The  values  of  the  parameters  A,  B,  C,  A,  r„,  r,,.  and  R,, 
are  those  of  Ref.  1. 
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to  explain  easily  the  well-known  experimental  fact 
that  and  ■*£,  for  Mn**3<f®  are  different  lines 
which  according  to  B,  C,  A  theory  are  degenerate. 
In  other  words,  certain  lines  which  are  accidental¬ 
ly  degenerate  in  the  fi,C,  A  theory  are  resolved  in 
the  improved  theory. 

The  present  theory  is  still  approximate  since  no 
consideration  has  been  given  for  the  spin-orbit  in¬ 
teraction  and  for  something  equivalent  to  what  is 
known  as  the  Trees  correction""*®  which  deals 
with  the  orbit-orbit  interaction  and  because  it  in¬ 
cludes  interactions  only  from  a  limited  number  of 
configurations.  With  extra  efforts,  however,  it  is 
possible  to  make  such  corrections  to  the  results 
calculated  in  the  present  framework.  Nonetheless, 
the  improvements  obtained  in  the  present  study 
especially  in  terms  of  generality,  are  significant. 
It  is  hoped  that  the  present  treatment  will  be  ex- 
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tended  appropriately  to  various  other  case.s’'  where 
the  experimental  results  are  in  poor  agreement 
with  the  theory.  Also,  it  is  expected  that  this  work 
will  provide  the  necessary  stimulus  to  experi¬ 
mentalists  to  observe  more  refined  and  abundant 
lines  in  the  spectra  and  to  the  theorists  to  make 
the  first-principles  calculations  to  explain  the  ob¬ 
served  results. 
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matrix  of  dimension  (43  x  43)  and  the  experimental 
energy  values  [5]  to  deduce,  in  units  of  cm"',  k',  — 

=  11,176,  Vt-y4=  11,931,  I',  -  Ks  =  647,  = 

2469,  Vt  =  6026,  K,  =  2725,  =  1584,  1^,  =  2164, 

Vio  =  6093  (and  Dq  =  571)  which  compare  with  the 
respective  free-ion  values  calculated  using  dementi's 
wave-functions  [8],  in  cm"'  15,328,  15,328, 0,  3958, 
8806,  5378.  1979,  1979,  7664  (with  the  calculated 
free-ion  value  Vi  =  195,437 cm"').  Again,  there  are  no 
values  available  for  f'j’s  from  other  sources  for  further 
comparison;  our  Dq  value,  however,  lies  close  to  the 
Dq  value  for  Mn^*  in  other  crystals  ( 1  ].  We  have  also 
obtained  assignments  of  the  lines  which  do  not  com¬ 
pletely  agree  with  the  tentative  assignments  [9]  given 
by  Yun  et  oL  [5);  in  particular,  the  second,  sixth, 
eighth,  and  ninth  lines  come  out  to  be  *Tig(*G), 
^TigCD),  ^AigCG)  and  *Ttg{*F)  instead  of '‘7’„(''C), 
*TtgCF)  and  ^Tigi^F),  respectively. 

In  the  present  treatment  we  have  not  taken  into 
account  the  width  of  the  observed  lines  in  our  calcu¬ 
lations  which  are  expected  to  change  our  numerical 
values  of  l^,'s  slightly  and  bring  about  better  agreement 
with  the  lines  not  involved  in  the  inverse  eigenvalues 
problem.  Also,  we  have  appropriately  assumed  that  the 
possible  departures  from  the  cubic  symmetry  and  the 
spin-orbit  interactions  in  these  systems  are  negligible. 
Though  the  theory  requires  many  experimental  lines, 
it  is  important  since  it  gives  information  about  the 
various  Coulomb  and  exchange  integrals.  The  present 
treatment  is,  however,  approximate  in  the  sense  that 
only  a  restricted  manifold  of  configurations  have  been 
incorporated. 

The  departure  from  the  cubic  symmetry  at  the 
impurity  site  due  to  the  F-center  has  been  analyzed  by 
Sibley  etaL  [10].  They  observed  the  nonuniform  shift 
of  the  lines  by  as  much  as  1000  cm"'  on  comparing 
the  peak  positions  for  the  Mn’*  transitions  in  KMnFa 
and  MnF2  with  those  of  the  irradiated  doped  KMgF3 
and  MgFj.  While  the  observed  shifts  could  be  due  to 
the  effect  of  the  F-centers,  they  are  partly  due  to  the 
differences  in  the  crystal  structure  parameters  of  the 
systems  analyzed.  In  view  of  this  one  may  infer  that 
though  some  effects  on  the  numerical  values  of  the 
presently  deduced  Coulomb  and  exchange  integrals 
due  to  the  presence  of  F-centers  are  expected,  it  may 
not  be  very  drastic  in  changing  our  present  results. 

It  is  hoped  that  this  work  will  stimulate  exper¬ 
iments  to  observe  abundant  and  sharp  lines  also  in 
many  more  crystals  with  different  transition  metal 
ions;  it  will  then  be  possible  to  obtain  information  as 
to  the  Coulomb  and  exchange  interactions  (together 
with  refined  values  of  A)  in  those  systems  following 
the  procedure  adopted  here.  Next  step  in  line  is  to 
develop  the  first-principle  theory  to  explain  these 
interaction  constants.  For  this  purpose,  the  theoretical 
approaches  [11]  which  treat  the  impurities  in  solids 


appropriately  would  be  very  helpful.  This,  in  turn  will 
assist  in  elucidating  the  electronic  structure  of  impurities 
in  solids. 

In  summary,  we  have  used  a  refined  analysis  of  the 
transition  metal  ions  in  solids  which  not  only  reveals 
the  important  differences  frcni  the  conventional  three- 
parameter  theory  as  regards  the  assignments  of  the 
optical  lines  but  also  provide.;  estimates  of  the  Coulomb 
and  exchange  integrals  (along  with  the  refined  value  of 
A  or  \0Dq).  This  is  of  great  significance  since  such  an 
information,  to  the  knowledge  of  the  authors,  has  not 
been  available  from  any  other  source.  Also,  since  the 
present  treatment  is  based  entirely  on  the  group  theory 
in  conjunction  with  the  experimental  data  it  is  expected 
to  give  correct  assignments  to  the  experimental  spectral 
lines. 
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able  to  derive,  from  straightforward  but  tedious  group 
theoretical  considerations,  the  general  formula  which 
relates  a  diagonal  element  of  one  with  the  correspond¬ 
ing  one  of  the  other,  as 

<rr'"(5,  r,  r,)5rAf^r  I  •  i  r, 

X  (5jrj)5rAf,Afr>  -  <t2"'(S,r,)e"(Sjr2)5rAf.Afri  • 

•U2"'(5,r,)e"(S2r2)5r3/.Afr>  =  (3-m)f',+ 


being  fixed  to  identify  the  ground  level  with  the  zero 
energy);  the  remaining  line  can  be  used  for  checking  the 
validity  of  the  treatment.  For  the  inverse  eigenvalue 
problem  we  have  derived  the  r/ ’-electrostatic  matrices 
(appropriate  to  Co’*)  containing  the  unknown  param¬ 
eters  Vi.  Denoting  the  electrostatic  matrix  by  Af  and 
the  energy  values  by  h\  one  writes, 

MX  =  t'X 


+  4(3  -  m)V2  +  (24  -  4m  -  641) Vi  +  V4 

+  (6-3n)Fs -(-(5m-  (2m  +  3n  -  i2)V.,  + 

-t-(2m-l-3M-  12)^  Fg-2(3-m)F,o  (1) 

where  Ps  are  the  irreductible  representations  of  the 
cubic  group  and  Si's  are  total  spin  quantum  numbers. 

The  importance  of  this  formula  arises  from  the  fact 
that  all  diagonal  elements  of  a  complementary  configur¬ 
ation  relative  to  the  configuration  d"  are  no 
more  altered  by  the  same  amount,  as  it  is  in  case  of  the 
simplified  theory  assuming  “pure"  rf-orbitals  and  using 
the  Racah  parameters. 

Knowing  the  electrostatic  matrices  we  now 
encounter  the  problem  of  determing  F,*s.  If  sufficient 
number  of  experimental  lines  are  available,  it  is  possible, 
though  cumbersome,  to  determine  the  P,’s  by  solving 
the  relevant  inverse  eigenvalue  problem.  By  analysis 
one  recognizes  that  in  the  free-ion  limit,  the  parameters 

^21  y*  and  Vs  contain  the  Slater-Condon  factor 
(or  the  Racah  parameter  A ).  Thus,  it  would  be 
appropriate  to  use  one  of  the  four  parameters  (as  Fg  or 
A  in  the  free-ion  case)  to  fix  the  energy  scale  to  facilitate 
comparison  with  the  spectral  lines. 

Next,  we  turn  our  attention  to  the  interpretation 
of  the  optical  absorption  lines  as  observed  in  MgF2  .Co’* 
and  MgF2:  Mn’*  by  Yun  et  ai  [5].  For  the  related 
research  works  one  may  consult  also  the  review 
article  by  Ferguson  [ll.  In  case  of  MgF2:  Co’*  the 
spin-allowed  bands  are  those  of*Ajg,  *T2g  and  *Tig 
states  and  the  other  transitions  are  spin-forbidden. 

The  spin-forbidden  transitions  have  .been  observed,  as 
mentioned  earlier,  by  lifting  the  forbiddenness  through 
exchange  interaction  with  the  color-centers.  The 
observed  energy  values  and  their  tentative  assignments 
based  on  the  previous  simplified  theory  (which 
effectively  assumes  “pure"  </-type  orbitals)  have  been 
tabulated  in  Table  I . 

In  case  of  MgF2:  Co’*  there  are  elevent  well- 
defined  observed  lines  (including  *T\g{*F)  as  the 
ground  state)  out  of  which  ten  lines  are  used  for  the 
solution  of  the  inverse  eigenvalue  problem  for  deter¬ 
mining  the  ten  parameters  V^,  V3 . . .  Kig,  A  (with  F, 


where  X  stands  for  the  eigenfunctions.  For  the  inverse 
eigenvalue  problem  one  prescribes  as  many  number  of 
eigenvalues  as  the  number  of  unknown  parameters  and 
solves  for  the  parameters,  obtaining  in  the  process  also 
the  remaining  eigenvalues.  The  dimension  N  of  the 
matrix  varies  with  the  occupancy  of  the  t/-orbitals;  A'  is 
20  for  Co’*  which  is  a  d’  system  and  it  is  43  for  Mn 
Mn’*(z/*).  The  present  inverse  eigenvalue  problem  is 
indeed  very  difficult  to  solve  since  the  number  of  eigen¬ 
values  prescribed  is  less  than  the  dimension  of  the 
n  itrix.  However,  we  have  succeeded  in  generalizing 
the  method  given  by  Bohte  [7]  to  solve  the  problem  in 
our  case. 

The  solution  of  the  inverse  eigenvalue  problem 
yields  the  unknown  parameters  which  have  been 
listed  in  Table  2  for  MgF2:Co’*.  The  calculated  energy 
eigenvalues  and  the  assignments  have  been  given  in 
Table  1  and  compared  with  the  observed  values  and  the 
previous  tentative  assignments  of  Yun  et  al.  [5]  based 
on  the  “pure"  (f-orbital  theory.  The  perusal  of  Table  1 
reveals  the  important  differences  that  the  previous 
assignments  ’Tijl’C),  ’zl|g(’C)and  should 

respectively  be  ’r2g(’C),  ’r2g(^C)and  ’i4|g(’G).  Also 
the  first  ten  energy  values  agree  with  the  observed  values 
as  expected  and  the  remaining  energy  value  (’Fg(’//)) 
lies  very  close  to  the  observed  value.  Here  we  need  to 
emphasize  that  the  present  procedure  concerns  with 
solving  the  inverse  eigenvalue  problem  which  is  basically 
different  from  the  procedure  of  fitting  with  the  para¬ 
meters  usually  adopted  in  this  field.  In  Table  2  the 
Coulomb  and  exchange  integrals  calculated  for  the  free- 
Co’*  ion  using  dementi’s  SCF— HF  wavcfunctions  [8] 
have  also  been  listed  for  comparison.  We  have  used  V, 
for  fixing  the  lowest  state  at  zero  energy  and,  conse¬ 
quently,  Vi,  Vn  and  Vs  have  been  shown  in  Table  2 
with  respect  to  F, .  The  calculated  value  of  F|  for  the 
free-ion  case  is  217,375  cm”’ .  It  is  interesting  to  note 
that  the  derived  values  of  the  integrals  have  changed  as 
much  as  50%  from  the  free-ion  values.  For  further  com¬ 
parison  the  values  of  the  Coulomb  and  exchange  integrals 
Fj’s  are  not  available  from  any  other  source.  However, 
the  value  of  Dq  equal  to  —  785  cm”'  has  been  reported 
for  Co’*  in  fluorine  environment  (in  KC0F3)  which  lies 
close  to  our  value  of  —  725  cm''  (Table  2). 

As  for  MgF2:Mn’*  we  have  used  our  £/*-electrostatic 
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Table  I.  Tabulation  of  the  observed  [5|  and  calculated  energy  values  (by  solving  the  inverse-eigenvalue  problem  as 
described  in  the  text)  for  MgE^rCo^*.  The  previous  tentative  group  theoretical  assignments  of  the  levels  have  also 
been  listed  along  with  the  present  assignments 


Previous 

tentative 

assignments 

Observed 
energy  values* 

(cm"') 

Calculated 
energy  values® 

(cm"') 

Present 

assignments 

7,519 

7,518.9 

*T2g(‘*F) 

11,905 

1 1 ,904.9 

%eG) 

15,385 

12,205.1 

^T.gi^G) 

17,452 

15,385.0 

*A^g(*F) 

^T^gi^G) 

18.519 

17,452.1 

^T^gi^G) 

*T,g(*P) 

19,305 

20,921 

19305.0 

^T,g(*F) 

^A,g(^G) 

20,833 

20,832.9 

^T,g(Hi) 

^T,geH) 

21,739 

21,739.0 

^A.g^G) 

^T2g.  ^T,g(^H) 

23,095 

23,094.9 

^T^gi^H) 

23,529 

23,529.0 

^T,g(^H) 

%eH) 

26,968.9 

•  With  respect  to  *Tig(*F)  as  the  lowest  level. 


Table  1  List  of  the  Coulomb  and  exchange  interaction  constants  for  Co**  in  MgFj  obtained  by  the  procedure 
decribed  in  the  text.  The  corresponding  values  calculated  by  using  free-ion  SCF— HF  wavefunctions  have  also  been 
given  for  comparison 


Coulomb 

exchange 

and 

crystal  field 
parameters 

V\  -  f'j 

F.-F, 

Present 

values* 

(cm"') 

13,034 

13,728 

Free  ion 
values^ 

16,920 

16,920 

(cm"‘) 

K, 


-2433  2725  8491  3219 

0  4379  9724  5932 


Vs  V,  K,o  Dq 


1928  1259  6257  -  725 

2189  2189  8460 


*  From  the  solution  of  the  inverse-eigenvalue  problem  as  explained  in  the  text. 

**  Using  Qementi’s  wavefunctions  [8] ;  the  calculated  value  of  Vi  for  free-Co**  comes  out  to  be  2 1 7,375  cm"' . 


therefore,  the  derived  electrostatic  matrix  elements 
contain  10  independent  Coulomb  and  exchange 
integrals  [4]  instead  of  the  three  Coulomb  repulsion 
parameters  (A,B,  and  C).  The  ten  independent 
parameters  are:  <£»7l '  <^11  •  lc{>, 

<efi  •  lef),  <eei- 100),  <001  •  |ee>,  <00\ •  iw),  <fle| •  ftm), 
<6n\  ’  l|f>  and  <|{|  *  Irp?)  which  will  be  referred  to  as 
Ki,  Kj  . . .  and  Vio,  respectively.  The  Vfs  are  the  two- 
electron  matrix  elements  with  the  “dot”  representing 
the  interaction  operator  eV/’n,  rij  being  the  distance 
between  the  two-electrons.  The  0,  e,  q  and  j;  stand 
for  the  one-electron  orbitals  in  the  cubic  field 


representation.  They  are  not  the  “pure”  cf-orbitals  if  the 
transition  metal  ions  concerned  are  embedded  in  solids 
since  the  solid  state  effects  perturb  the  free-ion  d- 
orbitals.  Remaining  two-electron  integrals  which  involve 
0,  e,  I,  Tj  and  J:  functions  can  be  shown,  from  the  group 
theoretical  arguments,  to  be  related  to  the  K,-’s. 

Since  the  electrostatic  matrix  elements  are  large  in 
number  and  contain  many  terms,  we  do  not  intend  to 
list  them  here  in  terms  of  the  Kj’s.  The  nondiagonal 
matrix  elements  of  the  complementary  states 
come  out  to  be  the  same  as  of  the  states  t-^e".  How¬ 
ever,  the  diagonal  elements  are  different.  We  have  been'" 
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Making  use  of  Racah’s  irreducible  tensor  operator  technique  a  refined 
treatment  for  the  3<f-transition  metal  ions  in  crystals  has  been  presented. 

The  refined  theory  has  been  applied  to  analyse  the  recently  observed  fine 
optical  spectra  of  MgFj.Co^*  and  MgF2:Mn**.  Important  differences 
have  been  noticed  in  the  assignments  of  the  optical  lines  obtained  from 
the  present  analysis  and  from  the  conventional  three-parameter  theory. 

First  numerical  values  of  the  exchange  and  Coulomb  integrals  in  these 
systems  have  also  been  given. 

THE  ELECTRONIC  ABSORPTION  SPECTRA  of  the  In  spite  of  that,  the  B,  C,  A-theory  has  served  useful 

transition  metal  ions  in  crystals  have  been  the  subject  purpose  owing  to  the  following  two  main  reasons.  First, 

of  extensive  experimental  and  theoretical  studies  [1].  the  experimental  data  were  usually  scanty  with  the 

The  first  application  of  the  theory  developed  for  com-  result  that  a  simple  theory  involving  only  a  few  param- 
plex  ions  in  crystals  incorporating  effects  of  the  electro-  eters  was  good  enough  to  analyse  the  spectra.  Second, 
static  field  of  the  environment  on  the  free-ion  energy  no  refined  theoretical  expressions  were  available  for 

levels,  was  made  by  Finkelstein  and  Van  VIeck  [2]  in  better  interpretation  even  when  many  spectral  lines 

1940.  Since  then  many  researchers  [1 , 3-S]  have  were  provided  by  the  experiments.  Recently,  Sibley 

extended  our  knowledge  of  the  electronic  structure  of  and  co-workers  [5]  and  others  have  performed  excel- 

the  transition  metal  ions  in  solids  considerably.  lent  experiments  in  which  they  have  been  able  to 

The  interpretation  of  the  experimental  spectra  for  observe  many  sharp  lines  of  the  transition  metal  ions  in 
the  iron-group  ions  in  crystals  has  hitherto  been  con-  crystals  by  lifting  the  spin-forbiddenness  of  the  exci- 
fined  to  the  use  of  the  Kacah’s  Coulomb  repulsion  tations  through  exchange  interactions  with  color  centers, 

parameters  [6]  B  and  C  and  the  crystal  field  splitting  These  experiments  have  furnished  not  only  the  spin- 
parameter  A.  The  parameters  [6]  B  and  C  were  origin-  quartet  states  but  also  many  hitherto  unknown  spin- 

ally  involved  in  the  theory  developed  by  Racah  to  doublets.  From  these  results  we  have  been  encouraged 

explain  the  spectra  of  free-ions,  in  particular,  the  free-  to  promptly  reOne  the  method  and  use  it  to  analyse  the 
iron-group  ions.  Later  on,  the  theory  was  extended  to  experiments.  In  the  following  we  first  outline  the 
explain  the  optical  spectra  of  the  transition  metal  ions  method  briefly  and  apply  it  to  MgF] :  Co^^  and 

in  solids  incorporating  the  effects  of  the  electrostatic  MgFi:  Mn^*  where  abundant  experimental  data  [S] 

crystal  fields  on  the  ions  assuming  that  the  outermost  have  recently  been  obtained, 
unfilled  electrons  were  “purely”  of  (/-character.  The  present  theoretical  treatment  is  well-known 

specially  in  their  angular  behavior,  as  in  the  case  of  and  is  based  on  the  strong  field  scheme  since  the  ions  of 

free-ions.  Though  this  assumption  simplifies  the  theory  the  Hrst  transition  series  in  many  of  their  compounds 

considerably  and  makes  the  interpretation  very  con-  belong  to  this  scheme.  The  one-electron  orbitals  are  now 

venient,  it  is  basically  faulty  since  the  wavefunctions  sets  of  orbitals  which  form  bases  for  irreducible  repre- 

of  an  ion  undergo  changes  when  it  is  embedded  in  a  sentations  of  the  symmetry  group  of  the  crystal  poten- 
crystal.  The  inaccuracy  of  the  theory  is  also  revealed  tial. 

by  the  fact  that  it  does  not  work  well  particularly  Hie  matrices  of  the  electrostatic  interaction 

when  more  spectral  lines  are  involved  in  the  interpret-  between  the  electrons  have  been  determined  by  utilizing 

ation  [2].  Also,  in  these  circumstances,  not  only  that  Racah’s  irreducible  tensor  operator  formalism.  This 
the  fitting  parameters  whitdi  it  yields  are  rough  in  esti-  method  is  not  new  and  has  been  used  previously  by 

mates  but  they  (particularly  B  and  C)  lose  their  real  Tanabe  and  Sugano  [3]  for  deriving  the  electrostatic 

physical  significance.  Moreover,  the  poor  agreement  of  matrix  elements  in  terms  of  the  Racah  parameters  A ,  B 

the  calculated  energies  with  the  observed  spectral  lines  and  C  with  the  simplifying  assumption  that  the  angular 

casts  doubts  on  the  correctness  of  the  assignments  of  parts  of  the  one-electron  orbitals  are  of  "pure"  (/-type, 
the  lines  predicted  by  the  theory.  The  present  derivation  is  not  subjected  to  this  restriction; 


